home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / dpptrs.z / dpptrs
Text File  |  1998-10-30  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          DDDDPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPTRS - solve a system of linear equations A*X = B with a symmetric
  10.      positive definite matrix A in packed storage using the Cholesky
  11.      factorization A = U**T*U or A = L*L**T computed by DPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPPTRS( UPLO, N, NRHS, AP, B, LDB, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDB, N, NRHS
  19.  
  20.          DOUBLE         PRECISION AP( * ), B( LDB, * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DPPTRS solves a system of linear equations A*X = B with a symmetric
  24.      positive definite matrix A in packed storage using the Cholesky
  25.      factorization A = U**T*U or A = L*L**T computed by DPPTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangle of A is stored;
  31.              = 'L':  Lower triangle of A is stored.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      NRHS    (input) INTEGER
  37.              The number of right hand sides, i.e., the number of columns of
  38.              the matrix B.  NRHS >= 0.
  39.  
  40.      AP      (input) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  41.              The triangular factor U or L from the Cholesky factorization A =
  42.              U**T*U or A = L*L**T, packed columnwise in a linear array.  The
  43.              j-th column of U or L is stored in the array AP as follows:  if
  44.              UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j; if UPLO =
  45.              'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
  46.  
  47.      B       (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS)
  48.              On entry, the right hand side matrix B.  On exit, the solution
  49.              matrix X.
  50.  
  51.      LDB     (input) INTEGER
  52.              The leading dimension of the array B.  LDB >= max(1,N).
  53.  
  54.      INFO    (output) INTEGER
  55.              = 0:  successful exit
  56.              < 0:  if INFO = -i, the i-th argument had an illegal value
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.